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Abstract 

We investigate the origin of intermittency for multiparticle distribution in mo- 
mentum space, following the idea that there is a kind of power law distribution of 
the space-time region of hadron emission. Using the formalism of current ensam- 
ble model to describe boson sources we discuss intermittency exponents for the 
coherent and incoherent ( chaotic) particle production scheme. 
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1 Introduction 



Recently, several experiments found that the phenomenon of in- 

termittency is dominated by very short range correlations between 

momenta of identical hadrons. As it is well known, the HBT correlations 
reflect the size and shape of the space-time region from which the observed 
identical particles are emitted . Remembering that intermittency is equiv- 
alent to a power law dependence of hadronic correlation functions ( mul- 
tiplicity distributions ), one may conclude that the power law dependence 
must also be present in the distribution of space-time shapes and sizes of the 
region of hadron emission The first analysis of this problem was done 
in Ref. using the HBT formalism for purely incoherent hadron source. 
The considered description of hadron emission was simplified : totally in- 
coherent production from the space region was assumed, and therefore the 
single particle distribution could not be adequately described. Nevertheless, 
intermittency exponents obtained in this model linearly increased with the 
rank of multiparticle distributions. 

In order to correct description of the single particle distribution the pion 
source should not only have a power-law density profile but also some space- 
time correlations between particle emitting points. In this paper we would 
like to consider the model for HBT correlations worked out by Gyulassy, 



Kauffmann and Wilson in Ref. 10 1, known as the covariant current formal- 



ism or the current ensemble model. It has been often used to describe many 



different aspects of multiparticle production ( see eg. Ref. [T^, |]T3[, for review 
Ref. 0] and references therein ) . Covariant current formalism is a generaliza- 
tion of the HBT formalism formulated in the language of field theory. Pion 
emitting regions are described as "currents" with different phases, and incor- 
porated into the pion field equation in the form of current ensemble. Each 
current is localized around some point called a collision site, and the collision 
sites are distributed in space-time with some probability density. Pion emis- 
sion depends on current phases: if phases are identical, particles are emitted 
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coherently, if they change randomly, the pion production will be incoherent. 

Current ensemble model suits our purposes very well. On the one hand, 
it is a natural continuation of the treatment presented in The quantum- 
mechanical source density is replaced by the power-law distribution of col- 
lision sites. The case of coherent and incoherent emission can be described 
easily. On the other hand, current ensemble formalism operates already on 
pions ( hadrons). This fact removes one serious difficulty which plagues the 
interpretation of intermittency in terms of parton distribution. The problem 
is that, even if one explains the intermittent behaviour of the parton system, 
the experiment observes hadrons, not partons. Consequently, one should take 
into account the hadronization process. Usually, one invokes at this point 
the principle of parton-hadron duality, the phenomenon which is in fact not 
understood. However, if intermittency is caused by the HBT phenomenon, 
there is no problem of hadronization and the decay of resonances: the HBT 
effect works already on hadron level, reflecting the structure of the source 
density. And so does it in the current ensemble model. 

Alternative approaches to the calculation of Bose-Einstein correlations 
( for review see |^ and references therein ) incorporate also the reaction 
dynamics. Hence the current ensemble model can be considered as a plausible 
" first approximation" concentrated only on the effects of source geometry 
on the shape of correlation functions. 

In our previous studies of the current ensemble formalism |1T| we have 
already proved that assuming purely coherent production of particles, mul- 
tiparticle distributions are infiuenced by the scaling space-time structure of 
bosonic source, and give observed intermittent results if the source number 
is large enough. 

In the present paper we analyze the case of incoherent ( chaotic) parti- 
cle production. We apply the method of random phase averaging described 
in Ref. [0 and to generate incoherence. This formalism is presented 
in section 2. In section 3 we explain briefiy how to measure intermittency 
parameters. In section 4 we incorporate the scaling dependence of the pro- 
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duction region into the model. At first we discuss the simple power law 
distribution of collision sites following and prove that it does give ob- 
servable intermittent result. Later we show that the obtained results can be 
generalized for other source distributions showing the scaling behaviour. At 
the end of this section we describe the special case of the fixed source number 
for correlation functions. 

The aim of our discussion is to prove that in a well working physical 
model of boson production like the current ensemble model, multiparticle 
distributions in momentum space can be strongly influenced by the scal- 
ing behaviour of the space-time structure of bosonic source, and can show 
intermittent behaviour. 



2 Classical current formalism 
2.1 Pion fields 

In this section we would like to recall the main ideas of current ensemble 
model |T^. To obtain the final pion state produced by a classical current 
source one should solve the following field equation for the scalar pion field 
$(a;) with the source current operator J(x) : 

(9^a'^ + m^)$(x) = J(x) (1) 

In principle, the source current J{x) is an operator coupled to the pion fields, 
and treating it as a complex space-time function is only an approximation. 
Because we will not specify the conditions of production process it is difficult 
to define when this approximation could be used. We assume simply that 
we are allowed in our model to replace the pion current operator J{x) by its 
expectation value [p!0| . 

Solution of (1) gives the coherent final multipion state : 

I $ >= e-'^expii J £kJ(k)a^(k)) | >, < $ | $ >= 1 (2) 
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where a(k), a+(k) are creation and annihilation operators, n is the average 
pion multiphcity and J{k) is the on- mass-shell Fourier transform of J{x) [0], 
[ |T0|] . Pion density matrix constructed from (2) is : 

p=|$><<|.| (3) 

Multipion- inclusive distributions Pm.{ki, . . . , km) are then defined by the for- 
mula : 

1 d^"^a(n . . . tt) 

Pm(ki, . . . , km) = = Tr[p^a'^{ki) . . . a+(km)a(kni) • • . a(ki)] 

(4) 

2.2 Pion source as current ensemble 

Now we are ready to consider the following description of pion emission: 
pions are produced in space-time centers Xi, . . . ,X]\f, distributed with some 
probability density p{x). These centers can be for example an effect of N 
separate collisions "producing" pions [|l^. In this picture, the total pion 
source J{x) would thus be a sum of N different currents Ji{x) : 

N 

.Jix) = Y.Mx)e'^' (5) 

i=l 

where we allowed the possibility that the individual currents have different 
phases. It reduces easily to the purely coherent case described in |TT| assum- 
ing that : 

Lpi = . . . = LPN = f (6) 

We shall consider the situation when the Ji{x) depend only on distance from 
the individual collision site Xj. It means that, if a collision centered in x = 
is parametrized by j^ix) we get : 

N 

J{x) = Y.Ux-Xi)e'^^ (7) 

i=l 
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and the on-mass-shell Fourier transform is given by : 



N 



J(k) = i^(k) J2 e''^'exp{iujkti - ik^i), ujk = + ml (8) 



i=l 



2.3 Multipion distributions for current ensemble. The 
coherent and incoherent case. 

The density matrix p^^ averaged over the number of sources N and positions 
Xi of the pion sources is : 

p^ = J2 PiN) I d^xi . . . d^XMp{x{) . . . p{xn) I $ X $ I (9) 

where P{N) denotes the probability to find exactly N pion sources, and p{x) 
is the probability density to find the source placed at the point x. The density 
p{x) is normalized to 1. Using the formulae (2), (4) and (9) one can obtain 
the inclusive multipion distribution Pm(ki, • • • , km) in the form : 



P^(ki, . . . , k^) = ^ P{N) j d^x, . . . d^XNpixi) . . . p{xn) I J(ki) 1^ . . I J(k^) |2 

(10) 

Substituting J(k) in (10) in the form (8) one obtains : 

P^(ki, . . . , kn,) =1 j^(ki) P . . . I j^{k^) P ^(^) / ^"^^1 • • • d^XNp{xi) . . . p{xn) 

N •' 

N N 

il = l «2m=l 

One should remember that in the scalar product of 4-dimensional vectors ki 
and Xj the first component of kj is equal to as a result of the on-mass- 
shell Fourier transform taken in (8). For current phases satisfying the relation 
(6) we recover multiplicity distributions for the coherent particle production 
from [^. If the sources add incoherently, we should also average (11) over 
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the set of phases. This method of incoherence generation has led to the name 
"chaotic" production because we assume that the current phases ipi fluctuate 
randomly from to 2ti. After phase averaging we obtain the formulae for 
the single- and double pion distribution in the form : 

Fi(ki)=| J.(ki) |2<iV> (12) 

P2(ki, ks) =1 j.(ki) n j,(k2) \''{<N^> + < N{N - 1) >| p(ki - ka) H 

(13) 

One should notice that for = 1 one gets the coherent field results Ref. 
[ pJ]) . The produced particles are then uncorrelated. Hence we exclude the 
case = 1 from our further investigations. 

3 Intermittency exponents 

Intermittency is equivalent to the power law dependence of hadronic moments 
in momentum space. The hadronic moments we would like to consider in this 
paper are factorial moments obtained by integration of multiparticle distri- 
butions and cumulants obtained by integration of correlation functions in the 
finite region of size 5. Generally, one studies the behaviour of the factorial 
moments which tests the whole n-particle distribution function. However, 
the cumulants have the advantage of testing the genuine n-particle correla- 
tions, and so it is always interesting to investigate their contribution to the 
higher order correlations (see e.g. fl^). 

Scaling dependence of the moments is described by the intermittency 
exponents. After integration of pion multiphcities Fm(ki, . . . , k^) in the 
finite region of size 6, one obtains a kind of series in 6 : 

Fm{S) = /'p^(ki,...,k^)rf3A;i...d3A;^ = ^a,(L5)-"^- (14) 
Jo J 

where aj are intermittency indices for the factorial moment Fm, cij are the 
weights and the length L is introduced for dimensional reasons. In experiment 
we observe the term dominating in (14). 
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So, for our purposes we deRne the intermittency exponent to be equal 
to Oj taken from the term dominating in (14). We should notice it need not 
to be the term with max/min(aj) because of the weights aj. In the similar 
way one can get the cumulant intermittency exponent Um- 

There is not much experimental data concerning the higher order mo- 
ments. The existing ones show intermittency exponents increasing with rank 
of multiparticle distributions [1^. The NA22 group measured also the third 
order cumulant with the preliminary result z/3 = 2^2 . 



4 Scaling behaviour of the distribution p{x) 

4.1 Scaling distribution of production centers 

In this section we show that the scaling distribution of collision sites p{x) 
( defined in (9) ) can produce the scahng of momentum distributions Pm- We 
consider the simplest example of the distribution p{x) with the assumption 
that all particles were produced at the same moment to : 

pix) = 5it-to)psi\^\) (15) 

and the space part follows a power law : 

P5(x) = I X |"~3^0 < a< 2,a ^ 1 (16) 

where the length L was introduced for dimensional reasons. 

The exact power law behaviour in (16) implies that the p{x) in (15) cannot 
be normalized. Following |^ we introduce a simple cut-off to get rid of the 
problem. It could be also done in a more elegant way but for our purposes 
of rather qualitative analysis it is enough to use the cut-off : 

P5(x) = I x 1"-=^ e(L- I x |)a(47r)-i (17) 

In this case L can be interpreted as a size of the pion production region. The 
assumption of simultaneous particle production in (15) allows to consider in 
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(11) and respectively in (12), (13) only 3-dimensional Fourier transforms of 
ps- Using the relations for 3-dimensional Fourier transform : 

P5(k) = / Sxe-'^^-psiyi) (18) 

P5(k = 0) = I d^xpsi-K) (19) 
one obtains the distribution in momentum space : 

ps(k) = a(L I k I)-" (^^^^ duu^'-'^smu (20) 
Jo 

One can observe the power law behaviour in (20) only for L | k |> 1. For 
L I k |< 1 the singularity will be cut off, and p(k) in (20) tends smoothly to 
J for L I k |— > 0. We notice also that P5(k) < 1 for any k. 

4.2 Leading terms for multiplicities in the chaotic source 
current ensemble 

Now we will discuss leading terms in (14) for pion multiplicities defined in 
(11). At first we consider the lowest order distributions (12) and (13). After 
integrating the one-, two-, and three-pion distribution from (15) in the region 
(0, 5) we get respectively : 

Fi{S) =<N> cn{S) (21) 

F2{6) =< > C2i(5)+ < A^(A^ - 1) > C22(<5) (22) 

F^id) =< > C3i{S)+ < N\N-1) > C32(^)+ < 7V(A^-l)(iV-2) > 033(6) 

(23) 

where arc results of integrating Pi(ki), P2(ki, k2), -P3(ki, k2, ks) obtained 
from (11): they depend on 6 but they arc independent of N. It is also impor- 
tant to notice that both terms in (22) which survive after phase averaging 
have the same rank of N. From (21) we obtain obviously ai = for Pi(ki). 
From the distribution P2(ki,k2) we will get the power singularity only 
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in C22{S). There are two singularities: C32 ~ and C33 ~ in F3 of 

the same rank in N, so the singularity 5^^" will dominate here. Hence the 
intermittency exponent is : 

/2 = /s = 2a (24) 

Now we will analyze the general form of multipion distribution Pm(ki, • ■ • , km) 
to get the dominating term there. At first let us consider terms which can 
survive in (11) after phase averaging. The answer will be obtained from the 
analysis of the expression : 

where the average is taken over the set of phases The exponent product 
in brackets should be equal to 1. All the surviving terms have the same rank 
A^™ ofN. 

So we are ready to consider leading terms in Pm(ki, ■ ■ ■ ,^m) for any 
m > 1. Because all the terms have the same rank of N, and ps(k) shows the 
power law behaviour only for L6 > 1 a.s mentioned in the previous section 
we conclude that the term with the smallest intermittency exponent will 
dominate and therefore : 

fm = 2a (26) 

One can check this result holds on for both large and small N. It means 
that intermittency exponents for factorial moments will not change with the 
rank of multiplicity for the chaotic production. The result for scaled factorial 
moments will be the same. 

The expression we have got above differs from intermittency exponents 



obtained in [|TT| for the coherent source ensemble. In |[TT| the intermittency 
exponents fulfilled the relation : 

fL = 2ma (27) 

The calculations we have done to get requested the assumption N to be 
large enough. For small N the formulae were very complicated, and it was 
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difficult to see how they could provide the simple power law behaviour as 
observed in experiment. There were also no scaling in the scaled factorial 
moments. 

The result (26) can be actually derived under much less restrictive condi- 
tions. To see this, let us ffist formulate general conditions that the function 
p{x) must fulfill to be considered as a probability distribution and to show 
the scaling behaviour : 

P5(x) > (28) 

psix)d^x = 1 (29) 

P5(k) oc const{L \ k l)^" in someintervalL \ k |g {a,b) (30) 

One can notice that conditions (28), (29) give the following inequality for the 
3-dimensional Fourier transform : 

Ps{k) < ps{k = 0) = 1 for any k (31) 

So we can generalize the results obtained for the collision sites distribution 
p{x) defined in (15), (16). For any function p{x) which fulffils conditions 
(28), (29), (30) intermittency exponents behave like (26). 

4.3 Correlation functions 



Analyzing the results from Ref. |]TT| for the coherent source ensemble one can 
observe that intermittency exponents (27) calculated from multiplicities are 
equal to intermittency exponents obtained from correlation functions. For 
the chaotic production we can make an interesting remark. If the source 
number N is fixed - in our formalism it means that N is given by a "spike" 
distribution P{N) and < >~ A^, the scaled correlation functions look 
like : 

C2(ki, ka) = I /9(ki - ka) P (32) 

C3(kl,k2,k3) = ^^^^^^^^^^[p(ki-k2)p(k3-ki)p(k3-k2)+C.C.] 

(33) 
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It seems to be possible to generalize easily the above results for correlation 
functions of any rank m > 1. Then intermittency exponents will grow with 
the rank m following the rule : 

I'm = ma (34) 

identically with the behaviour of intermittency exponents calculated in Ref. 
[^] . And for multiplicities we get as usually intermittency exponents following 
(26). There is no contradiction here. The term with the exponent (34) 
appears also in the multiplicity but it will be small compared with the term 
with the smallest exponent (26). So correlation functions seem to be a better 
tool to investigate the scaling properties of the pion production. In the 
above case the correlation function simply extracts the term with the largest 
exponent. 

5 Conclusions 

We have investigated the relation between intermittency and the scaling be- 
haviour of the space-time distribution of collision sites in current ensemble 
model. Our conclusions can be summarized as follows : 

- there is a possibility of intermittency in current ensemble model provided 
that 

I. source number N is not fixed and : 

(a) one assumes a power law singularity in the distribution of collision sites 

(b) the particle production is totally coherent 

(c) the number of coherent sources N is large enough. Then intermittency 
exponents grow linearly with the increasing rank of multiplicities fol- 
lowing the Eq. (27). For small number of sources N intermittency is 
generally not observed. There are no well defined leading terms which 
can give intermittency exponents growing with the rank of moments/ 
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multiplicities. In the scaled factorial moments/cumulants scaling does 
not appear. 

II. source number N is fixed and : 

(a) one assumes a power law singularity in the distribution of collision sites 

(b) the particle production is incoherent 

Then the intermittency exponents for scaled and not scaled factorial mo- 
ments do not change with the rank of multiplicity. However, intermittency 
exponents for scaled cumulants follow the result obtained in Ref. ||^ for one 
incoherent source, i.e. the formula (34). 

- the above results only partially agree with the experiment. The inter- 
mittency exponents can grow with the rank of the multiplicity/correlation 
function, except the intermittency exponents for the factorial moments in 
the case of incoherent production. In this case correlation functions better 
detect scaling properties of the emission region than multiparticle densities. 
However, the result for cumulants = 2i^2 obtained by the NA22 group 
cannot be confirmed. 

I would like to thank Professor A. Bialas for many stimulating discussions, 
suggestions and a continuous interest in this work. I am indebted to Professor 
K. Fialkowski for his reading the manuscript and many helpful comments. 
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